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We study an exact rotating and expanding solution of the fluid dynamical model of heavy ion 
reactions, that take into account the rate of slowing down of the rotation due to the longitudinal 
and transverse expansion of the system. The parameters of the model are set on the basis of realistic 
3+1D fluid dynamical calculation at TeV energies, where the rotation is enhanced by the build up 
of the Kelvin Helmholtz Instability in the flow. 

PACS numbers: 25.75.-q, 24.70.+S, 47.32.Ef 


I. INTRODUCTION 

In peripheral heavy ion collisions the participant sys¬ 
tem has large angular momentum. In hydrodynamical 
model calculations it has been shown that this shear 
leads to a significant vorticity m- When Quark-Gluon 
Plasma (QGP) is formed with low viscosity [2H5], new 
phenomena may occur like rotation =Bj, or turbulence, 
which shows up in form of a starting Kelvin-Helmholtz 
Instability (KHI) [7f[9|. It was also shown 10] that in 
peripheral collisions even if the shear flow is neglected 
and the same boost invariant longitudinal velocity is as¬ 
sumed at all transverse points, in the Color Glas Conden¬ 
sate (CGC) model initial transverse flow develops and it 
contributes to angular momentum in the same direction 
as the angular momentum arising from the target and 
projectile motion in the participant system. 

Rotation in heavy ion collision has only recently been 
considered, our aim is to detect it using two particle cor¬ 
relation. 

The Differential Hanbury Brown and Twiss (DHBT) 
method has been introduced earlier in El- Previously 
the method has been applied to a high resolution Particle 
in Cell Relativistic (PICR) fluid dynamical model jl2] . 
Here we will look at how the DHBT can be used for the 
exact hydro model mm- Rotation in exact models 
have been investigated as well in m- 

We look at the values from the exact hydro model and 
determine the effect rotation has on the correlation func¬ 
tions (CF) for detectors at different positions. 


II. CORRELATION FUNCTION FOR EXACT 
HYDRO MODEL. 

The two particle correlation function for this model 
is found with the same method used in m where the 
source function is 
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difference is q = ( P 2 — Pi), u p is the 4-vector velocity of 
the source, is the normal of the freeze out hypersurface 
and T(x) is the temperature distribution. The density for 
the exact hydro model m is given by 

N(r p , r y ) = N b -^t exp(-r 2 /2R 2 ) exp(-r 2 /2F 2 ) (2) 

or using the scaling variables in the out (p, R), side 
(<p, 0), long ( y , Y) directions 

N{s p ,s y ) = N B ^exp{-s p /2)exp{-s y /2), (3) 

where s p = r 2 /R 2 and s y = r 2 /Y 2 . 

We use the finite size cylindrical shape source as de¬ 
scribed in eq. (10) of [14] 
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and the integral J(k,q) for this model will be 
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where k 0 = ^+ k 2 and q 0 = 


k ■ q 


w s is a weight 

function, w s oc k^a pi and the temperature profile is flat 
with a value of 250 MeV. 

We have the single particle distribution integral 
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and the correlation function is given by El 


where C n is a constant, is the average 4 vector mo¬ 
mentum of two pions, k = {pi +p 2)/2 and the momentum 


C(k,q) = 1 + 
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so the constants outside the integrals will cancel. 

The velocity consists of a radial expansion in the out 
direction, an angular velocity, where the rotation is in 
the reaction plane and also an axis-directed expansion in 
the side direction. The velocity in the out (p), side (ip), 
long-direction (y) is then 

v s — Spi Ruj yf sp , Y yj , ( 8 ) 

or in the x, y , z directions where x is the direction of 
the impact parameter, y is the longitudinal (rotation) 


axis and z is the (beam) collision axis. 

v s = [RJs~p sin(</j) + RujJlTp cos(i p), 

. . , ( 9 ) 

Yy/Sy, Ry/Sp cos (ip) - RujyfSp sin(^) J . 

The average transverse radius is R = VXZ, and we 
use this value when the exact model is studied. Using 
the values from Table [l] we show the correlation function 
as function of q = q out in Figures |T] a and b for two 
different configurations. 
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FIG. 1. (color online) Correlation Function C{k,q ou t) for the exact hydro model as function of q = q ou t, with 

(a) R = 2.500 fm, R = 0.250 c, Y = 4.000 fm, Y = 0.300 fm, lo = 0.150 c/fm,at t = 0.0 fm/c. (left figure) 

(b) R = 3.970 fm, R, = 0.646 c, Y = 5.258 fm, Y = 0.503 fm, oj = 0.059 c/fm, at t = 3.0 fm/c. (middle figure) 

(c) R = 7.629 fm, R = 0.779 c, Y = 8.049 fm, Y = 0.591 fm, u = 0.016 c/fm, at t = 8.0 fm/c. (right figure) 

The solid black line is for k = 0.2 fm -1 , the dashed red line is for k = 1 /m _1 and the dotted blue line is for k = 5 /m -1 


As the system size increases with time we get a narrower 
distribution in q out for the correlation function. We also 
notice a wider distribution for increasing values of k. 

III. RESULTS - DIFFERENTIAL HBT FOR 
EXACT HYDRO MODEL. 

Let us now ’’Event by Event” evaluate two correla¬ 
tion functions at two different k-vectors in the plane of 
rotation. The differential correlation function Ell is ob¬ 
tained by taking the difference between the correlation 
functions with e.g. at detector positions k + = (a, 0, b)k 
and subtracting the CF for a detector at k~ = (a, 0, — b)k 
(a 2 + b 2 = 1) in x , y, z coordinates, 

A C(k,q) = C(k + ,q)-C(k~ iq ) . (10) 

The integrals for this model cannot be given in ana¬ 
lytic form, so the CF and DCF need to be integrated 
numerically. The momentum difference vector q may 
point in different directions, and it is usual to use the 


qout, qiong, qside system of directions. In the present 
work we show the q out dependence of the Correlation 
functions only, where k || q out . 

Detectors 

As we can see from eq. (8) the rotation leads to an 
asymmetry, both the p and ip components of the veloc¬ 
ity v s depend linearly on yfsp, thus the flow velocities 
of the system at a given FO time have a velocity pro¬ 
file Vs^ = const. vi P K The corresponding momentum of 
the fluid motion has approximately the same character¬ 
istics, and this is indicated by the red full line in Fig. [5] 
This characteristic flow velocity distribution determines 
the final momentum distribution of the emitted particles 
where the contribution of a fluid element at a radial co¬ 
ordinate s p leads to a tangential momentum distribution 
centered around = Ruiy/sf,. The resulting momen¬ 
tum distribution of the emitted particles from different 
radii are sketched in Fig. [2j 

Although our spatial source configuration is az- 
imuthally symmetric, our phase-space configuration is 
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(fm/c) 

) (fm) 

(c) 

(c/fm) 

(fm) 

(c) 

(Rad) 

0 . 

4.000 

0.300 

0.150 

2.500 

0.250 

0.000 

3. 

5.258 

0.503 

0.059 

3.970 

0.646 

0.307 

8. 

8.049 

0.591 

0.016 

7.629 

0.779 

0.467 


TABLE 1. Time dependence of some characteristic param¬ 
eters of the fluid dynamical calculation presented in ref. E- 
R is the average transverse radius, Y is the longitudinal length 
of the participant system, ip is the angle of the rotation of the 
interior region of the system, around the y-axis, measured 
from the horizontal, beam ( z ) direction in the reaction, \x, z] 
plane, R, Y are the speeds of expansion in transverse and 
longitudinal directions, and oj is the angular velocity of the 
internal region of the matter during the collision. 


ferent values of the correlation function, and ideally we 
would place them along or near the direction of highest 
speed in the side direction which would usually be in the 
beam direction. 


f (Pp'P<p) 



not, because of the given direction of the rotation. As 
eq. ([5]) shows the CF depends on q ■ v s and q ■ x, thus 
reversing the direction of either q or v s will change the 
CF! For vanishing rotation the CF would be the same for 
different azimuthal angles (e.g. d z = ±6) if the spatial 
source is azimuthally symmetric. If there is no expansion 
the rotation does not change the CF integrals either. 

We use a detector placed at (x, y, z ) = (0.935,0,0.353), 
so that k = (0.935,0,0.353)fc and q = q out = 
(0.935, 0, 0.353)q out , so that both k and q are parallel, 
are in the reaction plane and are orthogonal to the rota¬ 
tion axis, y. 

If we look at the azimuthal integrals for k ■ v s for the 
detectors at k = (ak,0,±bk) in eqs. , we have the 
integrals below, which will have a non-zero difference for 
a ^ 6, R ^ 0 and u/0, both for k and q out - 

f 2 ” ( f ■ 

/ exp I k^/~s~p ( a[R sin(tp) + Rui cos(tp)] 

+c[R cos(<p) — Ruj sin(</>)]^ cLp 
^ J exp [k^fsp ^a[Rsin(<p) + Rlo cos(tp)] 

—c[I?cos(</ 2 ) — ffixsin^)]) ) dtp . (11) 


This is so because, although the spatial distribution does 
not depend on ip, the momentum space distribution de¬ 
pends on v v . 

In the case of a realistic dynamical configuration, espe¬ 
cially in the initial configuration, the local rotation veloc¬ 
ity component from the shear flow is maximal at higher 
distances in the ±x directions and pointing towards the 
±z direction. 


FIG. 2. (color online) The schematic phase space distribu¬ 
tion of the rotating and expanding source in the momentum 
space. The momentum of the expansion increases with the 
radius just as the momentum arising from the rotation. So, 
higher radial flow momenta correspond to higher rotation mo¬ 
menta also, as given by eq.Q and indicated by the red line. 
At constant p p the distribution peaks at the momentum p^ in¬ 
dicated by the red line; the Jiittner distribution, eq.|8|, is not 
symmetric, it is elongated towards higher momenta, see Sec. 
2.4.2 of [16]- The resulting thermally smeared distribution is 
indicated by the contour lines. 


In the exact model discussed here we can see that at 
small values of k there is little to no difference in the 
DCF, but the difference will increase for higher values. 
As the system grows in size the distribution becomes 
more narrow for the CF and the DCF will become smaller 
for larger values of the relative momentum q. 

For the initial time, t=0, the DCF is small and pos¬ 
itive, Figure [3] a, but for the later times the amplitude 
is larger and negative. Comparing Figures [3] b and c we 
see that as expected the peaks of the DCF are more to 
the left for larger size because R ~ 1/q at half width of 
the correlation function. We also see that the amplitude 
for DCF for Figures [3] b and c are about the same. As 
the system continues to increase in size we would see a 
decrease in the amplitude because of the lower 00 value. 
For higher values of the temperature we get a smaller 
amplitude in the DCF, or for smaller values of the tem¬ 
perature we get a higher amplitude. At late times, t = 8 
fm/c and larger size, the CF is much more narrow as 
shown in Figure [I] c. Notice that at low wave numbers, k 
the CF has several zero points. This is also reflected in 
the DCF at the same k value, see Figure [3] c. 


For different detector positions we would expect dif- 
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FIG. 3. (color online) Differential Correlation Function for the exact hydro model as function of q = q ou t, with 

(a) R = 2.500 fm, R = 0.250 c, Y = 4.000 fm, Y = 0.300 fm, uj = 0.150 c/fm at t = 0.0 fm/c. (left figure) 

(b) R = 3.970 fm, R = 0.646 c, Y = 5.258 fm, Y = 0.503 fm, u > = 0.059 c/fm at t = 3.0 fm/c. (middle figure) 

(c) R = 7.629 fm, R = 0.779 c, Y = 8.049 fm, Y = 0.591 fm, w = 0.016 c/fm at t = 8.0 fm/c. (right figure) 

Where the solid black line is for k = 0.2 fm _1 , the dashed red line is for k = 1 /m _1 and the dotted blue line is for k = 5 fm -1 

In (a) the solid black line is close to the axis. 


The DCF is dependent on the positions of the detec¬ 
tors as demonstrated in Figure [3] This dependence can 
be used to maximize the amplitude of DCF in a given 
configuration, based on previous theoretical estimates. 
Measuring a single CF at different azimuthal angles in 
the plane of rotation does provide the same C(k,q )~s as 
our model is azimuthally symmetric! Thus, the difference 
is caused by the measuring two CFs in one event with the 
same source function homogeneity area, but two differ¬ 
ent k-vectors at different azimuths with respect to the 
source area. This is also indicated by eq. (27) of ITl] . 

( 2 ^. \ 

—jr-®- 1 factor appears, which leads to 
the sensitivity on vector k. 


IV. CONCLUSION 

The model calculations show that the Differential HBT 
method can give a measure of rotation in this exact hydro 
model. The differential correlation function is dependent 
on shape, temperature, radial velocity and angular ve¬ 
locity. Also the detector position is important. 

If we eliminate rotation or the radial expansion the 
DCF vanishes in the model. It also indicates that using 
the estimated small rotation velocities, ui = 0.01 — 0.15 
c/fm, we get a DCF value approaching 2-3 %. 
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